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STRAIN  ENERGY  RELEASE  RATE  ANALYSIS  OF  CYCLIC  DELAMINATION 
GROWTH  IN  COMPRESSIVELY  LOADED  LAMINATES 

John  D.  Whitcomb 
NASA  Langley  Research  Center 
Hampton,  Virginia  23665 

SUMMARY 

Delamlnatlon  growth  In  compresslvely  loaded  composite  laminates  was 
studied  analytically  and  experimentally.  The  configuration  used  In  the  study 
was  a  laminate  with  an  across-the-wldth  delamlnatlon.  An  approximate  super- 
position  stress  analysis  was  developed  to  quantify  the  effects  of  various  geo¬ 
metric,  material,  and  load  parameters  on  mode  I  and  mode  II  strain  energy 
release  rates  Gj.  and  Gjj,  respectively.  Calculated  values  of  Gj  and  Gj^  were 
then  compared  with  measured  cyclic  delamlnatlon  growth  rates  to  determine  the 
relative  Importance  of  Gj  and  High  growth  rates  were  observed  only  when 

Gj  was  large.  However,  slow  growth  was  observed  even  when  G^  was  negligibly 
small.  This  growth  apparently  was  due  to  a  large  value  of  G^ j £*  (f^  ” 

INTRODUCTION 

In  composite  structures  subjected  to  compression  loads,- delaminations  can 
cause  localized  buckling  (fig.  1).  High  Interlaminar  stresses  at  the  edges  of 
the  buckled  region  often  lead  to  cyclic  delamlnatlon  growth  (herein  referred  to 
as  instability-related  delamlnatlon  growth) . 

The  objective  of  this  paper  Is  to  Investigate  the  mechanism  of  Instability- 
related  delamlnatlon  growth.  Figure  1  shows  the  configuration  used  In  the 
study  a  laminate  with  a  "through-width”  delamlnatlon.  This  configuration 
was  selected  because  It  Is  perhaps  the  simplest  configuration  that  exhibits 
Instability— related  delamlnatlon  growth.  Goals  of  the  Investigation  were: 

(1)  to  develop  and  use  an  approximate  superposition  stress  analysis  to  explain 
how  various  geometric,  material,  and  load  parameters  affect  Interlaminar 
stresses,  (2)  to  determine  the  delamlnatlon  growth  behaviors  predicted  by 
several  different  criteria  based  on  strain  energy  release  rates,  and  (3)  to 
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compare  analytical  calculations  with  experimental  observations  to  determine 
the  applicability  of  each  growth  criterion. 

Because  of  the  stress  singularity  at  the  end  of  the  delamination  (crack 
tip),  calculated  stresses  there  have  little  meaning.  Strain  energy  release 
rates  are  finite  parameters  which  characterize  the  intensity  of  the  stresses 
near  the  crack  tip.  Consequently,  in  the  following  discussion  strain-energy 
release  rates  will  be  used  to  characterize  the  severity  of  the  interlaminar 
stresses. 


a 

a 

Aa 

b 

C,  Cl,  C2,  C3, 
n,  nl,  n2,  Z 

D 


NOMENCLATURE 

half-length  of  delauiination  before  loading 
half-length  of  delamination  after  loading 
virtual  crack  closure  distance  used  in  strain  energy 
release  rate  calculations 
specimen  width 

arbitrary  constants 

bending  stiffness  of  the  buckled  region  given  by 


D 


k«l  ^ 


where  p  «  number  of  plies 

dj^,  dy  unit  load  solutions  for  displacements  near  crack  tip 

E^  Young's  modulus  for  ply  k 

Ell,  ^22*  ^33  Yo'ing's  moduli  of  unidirectional  ply.  The  subscripts  1, 

2,  and  3  refer  to  the  longitudinal,  transverse,  and 
thickness  directions  respectively. 


<o 


ORIGINAl.  FAG?-  F5 
OF  POOR  QUALIff 


I 


®12*  ®13»  ^23 
M 

N 

^A*  ^C’ 

Sa»  Sg,  Sf,,  Sjj 


t 

X,  y 
“k-l 

6 

/» 

6 

Vi2.  Vj_3,  V23 


unit  load  solutions  for  forces  at  crack  tip 
Mode  I  strain-energy  release  rate 
Mode  II  strain-energy  release  rate 

maximum  possible  value  of  Gj  for  current  delamination 
length 

shear  moduli  of  unidirectional  ply 
moment 

number  of  applied  load  cycles 

axial  loads  in  regions  A,  C,  aud  D  respectively 
remote  applied  compressive  load 

axial  stiffness  of  regions  A,  B,  C  and  D  given  by 
P 

S  -  b  X)  -  cl  ) 

k«l  K 

where  p  •  number  of  plies 
thickness  of  buckled  region 
rectangular  Cartesian  coordinates 
distance  from  top  surface  of  laminate  to  ply  *'k”; 
top  ply  is  ply  1 

lateral  deflection  at  x  *  -a  due  to  applied  load 
value  of  6  corresponding  to 
initial  lateral  deflection  at  x  *  -a 
Poisson's  ratios  for  unidirectional  ply 


ANALYSIS 


The  configuration  shown  in  figure  1  was  Idealized  as  a  two  dimensional 
plane  strain  problem.  Linear  and  nonlinear  finite  element  analyses  and  an 
approximate  superposition  analysis  were  used  to  calculate  strain  energy 
release  rates  for  the  two  dimensional  idealization*  The  nonlinear  analysis 
was  used  to  provide  reference  solutions  for  evaluation  of  the  approximate 
superposition  analysis.  The  linear  analysis  was  used  to  calculate  several 
constants  used  in  the  approximate  superposition  analysis*  The  nonlinear 
analysis  is  described  in  reference  2,  and  the  linear  analysis  is  simply  a 
linear  version  of  this  analysis. 

The  approximate  superposition  analysis »  the  procedure  for  calculating 
strain  energy  release  rates,  the  finite  element  models,  and  material 
properties  are  discussed  in  the  following  sections. 

Approximate  Superposition  Analysis 

Superposition  techniques  have  been  widely  used  in  linear  stress  analysis 
to  represent  a  complicated  problem  as  a  combination  of  several  simpler  problems 
x\pplication  of  the  principle  of  superposition  to  nonlinear  problems  first 
requires  a  transformation  that  results  in  a  linear  system* 

The  key  to  the  transformation  is  replacement  of  the  source  of  nonlinearity 
with  equivalent  loads  (fig.  2a  and  b) .  Because  of  syranetry  only  half  of  the 
configuration  is  considered*  The  buckled  region  (which  responds  nonllnearly 
due  to  significant  rotations)  is  replaced  by  the  loads  Pp  and  M,  the  axial 
load  and  moment  respectively  in  the  column  where  it  is  cut  (fig.  2b).  The 
new  configuration  is  linear,  with  three  nonlinearly  related  applied  loads 
Pt»  superposition  the  number  of  loads  can  be  reduced  to  two,  as 

illustrated  in  figures  2c  -  2e. 
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The  load  system  in  figure  2c  is  divided  into  the  two  load  systems  shown 
in  figures  2d  and  2e«  Because  is  calculated  using  rule  of  mixtures « 

the  load  system  in  figure  2e  causes  a  uniform  axial  strain  state  and  no  inter* 
laminar  stresses.  Consequently^  in  terms  of  interlaminar  stresses,  only  the 
load  system  in  figure  2d  (ie.  and  M)  need  be  considered*  Accordingly, 

in  the  current  study  involving  strain  energy  release  rates,  figure  2d  is 
the  linearized  equivalent  of  the  nonlinear  problem  in  figure  2a* 

The  appendix  describes  a  strength  of  materials  analysis  for  calculating 

equations  from  the  appendix  are 


(6  ^  266^)  D(S^  +  Sjj) 


16 


A  D 


6+6 


(1) 


P 


C 


M 


6 

2a^ 


n^D  6 


(2) 

(3) 


To  use  the  loads  (P^  -  and  M  in  a  two  dimensional  analysis  requires  that 
they  be  expressed  as  an  equivalent  distribution  of  tractions.  To  calculate 
this  distribution,  the  axial  strains  were  assumed  to  vary  linearly  through  the 
thickness  where  the  tractions  are  applied  (ie.  at  the  cut).  Intuitively,  this 
seems  to  be  reasonable  if  region  D  (fig.  2)  is  not  cut  too  close  to  the 
crack  tip.  The  validity  of  the  assumed  linear  variation  will  be  checked  later 
in  this  paper. 

Linear  finite  element  analysis  was  used  to  calculate  the  response  of  the 

linearized  configuration  in  figure  2d  to  unit  values  of  (P  ~  P  )  and  M. 

CD 

Because  the  configuration  is  linear,  the  solution  for  any  arbitrary  combination 
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■  ^D>  M  Is  simply  a  linear  combination  of  the  unit  load 

responses.  If  region  B  (fig.  2)  is  much  thicker  than  region  C,  the  unit 
load  solutions  are  very  insensitive  to  delamlnatlon  length.  In  the  current 
study  the  ratio  of  thicknesses  was  61  to  3.  Hence,  the  unit  load  solutions 
for  2a  -  25  mm  were  used  for  analysing  all  delamlnatlon  lengths.  Also  initial 
wavlness  of  the  buckled  region  does  not  enter  into  the  finite  element  analysis. 
Delamlnatlon  length  and  initial  wavlness  were  both  accounted  for  In  the 
strength  of  materials  analysis  In  calculating  (Pg  -  P^)  and  M,  equations 

(2)  and  (3)  respectively.  This  procedure  will  be  discussed  further  In  the 
next  section. 

Strain  Energy  Release  Rate 

The  virtual  crack  closure  method  (ref.  1)  was  used  to  calculate  mode  I 
and  mode  II  strain-energy  release  rates,  and  respectively.  The  forces 
transmitted  through  the  node  at  the  crack  tip  and  the  relative  displacements 
of  the  two  nodes  on  the  crack  boundary  closest  to  the  crack  tip  were  used  In 

the  calculation.  Equations  4  show  how  this  teclmlque  is  used  for  the  super¬ 
position  stress  analysis. 

<=:  ■  2^  [(fc  -  ]  [<Pc  -  P„)<iJ  .  wjj 

if  1  r  1 

In  these  equations  F^,  F^,  d^,  and  d^  are  the  unit  load  values  of  the  nodal 
forces  and  the  corresponding  relative  nodal  displacements  in  the  x  and  y 
directions.  (The  coordinate  system  Is  defined  in  fig.  2.)  The  superscripts 
1  and  2  on  the  unit  load  parameters  Identify  parameters  associated  with 
«nd  M,  respectively. 
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If  the  distance  is  small  between  the  crack  tip  and  the  nodes  used  to 
calculate  relative  displacements,  then  Fy/dy  -  Fy/d^  and  F^/d^  - 
Using  these  relationships  in  eqns.  (4)  results  in 

y  L  J 

1  <5) 

'II-2S7  [cc- 

In  the  results  and  discussion  section  it  %d.ll  be  shown  that  for  high  loads 
or  long  delamlnatlon  lengths,  Gj  is  aero,  i.e.,  the  crack  tip  closes  in  the 
normal  direction.  To  prevent  the  crack  faces  from  overlapping  (analytically) 
requires  the  addition  of  multipoint  constraints  on  the  crack  face  nodes.  Con¬ 
ceptually,  the  crack  face  nodes  are  connected  in  the  direction  normal  to  the 
crack  face  by  infinitesimal  springs.  These  springs  have  infinite  stiffness 
in  compression  and  zero  stiffness  in  tension.  To  determine  whether  to  select  zero 
or  infinite  stiffness  requires  solution  of  a  non  tear  contact  problem.  To 
include  the  contact  problem  directly  in  the  superposition  analysis  would 
severely  complicate  the  otherwise  simple  equations.  Therefore,  use  of  a  non- 
contact  analysis  to  approximate  Gjj  was  investigated. 

A  laminate  with  2a  -  76.2  mm  was  analyzed  using  two  different  approaches. 
First  contact  forces  were  ignored  (i.e.,  overlap  of  crack  faces  was  allowed). 

Gj  and  Gjj  were  calculated  using  equation  (5).  In  the  second  approach,  overlap 
of  the  crack  faces  was  prevented,  which  is  more  realistic.  Gjj  was  calculated 
using  equation  (5).  (Note  that  eq.  (5)  yellds  Gj  -  0  when  overlap  is  prevented.) 
Applied  loads  (P^)  ranged  from  14.8  kN,  which  corresponds  approximately  to 
initial  crack  tip  closure,  to  55.2  kN. 
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VJhen  crack  face  overlap  was  prevanted.  a  larger  value  of  G^  was 
calculated  than  when  overlap  va.  aUowed.  The  difference  i»  the  G^i  values 
..c....ed  vlth  W.  But  lu  11  th.  dlft.r.~.  «~  .PPro.i».t.l,  .,u.l 

t,  Oi  clculutud  uslug  th.  .pptulh  I.lch  11»»d  ctuclc  f.c.  puetup.  Bor 
e««ple,  foe  P,  -  55.2  M  the  co»t«:t  eoeljel.  yielded  G^j  -  *13  J/»  •  ““h 

ctech  fete  oyetlep  11-*.  <>1  «U 

The  .u.  of  theee  .1...  f  HI.  epptol-tel,  1.5  petcent  of  the  ~re 
ee.ll.tlc  .Olutlc  l.e..  G„  -  *13  3/.^ 

loece.  do  not  .lB.lflc«.tlj  Itee  the  totl  eteel.  -.etc  eeleM.  tele.  Hence. 
When  there  is  crack  tip  closure,  the  total  strain  energy  release  rate  frost  the 
non-contact  «,.l,.l.  ten  he  need  to  .Ppeo1«t.  G„  (Hch  1.  the.  th.  totl 

etteln  eneeg,  e.1.1.  ente.  .me.  G,  1.  Id«.tlc11y  teto). 

Finite  Element  Model 

A  typical  finite  element  nenh  toe  the  nollneee  eneiyele  Is  sHoun  m 
flguee  3.  B.C.U..  of  aymmetey  oly  hlf  of  th.  llente  .»  ~dole*- 

„eh  contemn  ei3  node,  md  7*0  toue-nod.  leopeemetelc  elenente.  Beduced 
integeetlon  ...  uaed  to  Impeov.  th.  peefoe-.nco  of  the  elm-nt.  In  «.d.lln, 
bending  defoematlon..  Been.,  th.  eot.tl«..  .e.  «lll  ox~Pf  *»  ■>«' 
huchled  eeglon.  the  nonllne.e  .te.m-dmpl.c»».t  eeUtlon.  «e.  need  only  foe 
th.  eeglon  y  >  0.  -.  <  »  5  -0.56  «.  *•  .ho».  m  flgne.  «b).  th. 
im..elmd  conflgueetlon  1.  th.  1  the  nonimce  conflgue.tlon  encept  th.t 

«,.t  of  the  huchled  eeglon  1.  e.«.y.d.  Accoedlngly,  th.  .-.h  need  In  th. 

im..e  «..iy.i.  »•«  Ef" 

deleted  part  of  the  buckled  region. 
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Materials  Properties 

The  material  studied  was  KARMCO  T300/5203*  graphlte/epoxy.  The 
unidirectional  ply  properties  were  assumed  to  be 

Ell  “ 

^22  "  ^33  **  GPa 
''X2  - ''13  -  '^23  •  0*21 
°12  "  ^13  "  °23  “  5.9  CPa 

Plane  strain  (i.e.,  -  0)  and  -  0  were  Imposed  to  calculate  the  2D 

properties.  In  regions  where  coarse  finite  elements  spanned  several  plies, 
laminate  theory  was  used  to  obtain  average  properties. 


♦Use  of  trade  names  or  manufacturers  does  not  constitute  an  official 
endorsement,  either  expressed  or  implied,  by  the  National  Aeronautics  and 
Space  Administration. 
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EXPERIMENTAL  PROCEDURE 

The  specimens  used  for  this  study  vere  fabricated  and  tested  by  Northrop 
Corporation.  (Details  will  appear  in  NASA  CR-166046,  "Performance  of  a  Quan¬ 
titative  Study  of  Instability-Related  Delaminatlon  Growth,"  by  R.  L.  Ramkunar.) 

A  cursory  description  of  the  experimental  procedure  is  given  herein. 

The  specimen  consisted  of -64  plies  of  T300/5208.  The  fiber  orientation  and 
stacking  sequence  were  [0^/(0/45/90/^5}^]^.  The  laminate  width,  b,  was  25.4  nt 
To  simulate  a  delaminatlon,  k^ton  film  was  uaed  to  prevent  bonding  over  a 
19  mm  length  between  the  third  and  fourth  plies.  The  ply  thickness  was 
assumed  to  be  0.14  mra.  Six  specimens  were  tested  in  fatigue  under 
compressive  constant-amplitude  loads.  Minimum  compressive  load  was  ten  per¬ 
cent  of  the  maxlwm  con^resslye  lo^f  ^;  The  load  frequency  was  10  Hz.  Delami- 
natlon  lengths  .  were  measured  wl^>a  microscope. 

,1?,  results  AND  DISCUSSION 

First,  the  accuracy  of  the  approximate  superposition  analysis  will  be 
evaluated.  Then  the  effect  of  various  parameters  on  and  G^^  vlll  be  con¬ 
sidered.  Finally,  the  experimental  observations  will  be  compared  with  the 
analytical  results.'  '  ' 

Evaluation  of  Approximate  Superposition  Analysis 

The  approximate  superposition  analysis  was  evaluated  by  comparison  with 
results  from  a  geometrically  nonlinear  finite  element  analysis.  Recall  that  a 

<i 

major  assumption  in  the  approximate  analysis  was  that  the  strains  vary  linearly 
through  the  thickness  where  (P^  -  Pjj)  and  M  are  applied.  Fig.  4  shows  the 
axial  strain  variation  through  the  thickness  at  s  ■>  -0.2  and  -0.7  mm  obtained 
using  the  nonlinear  finite  element  analysis.  Along  the  line  x  ■>  -0.7  mm  the 
strains  vary  alioost  linearly  for  the  three  applied  loads.  However,  closer  to 

X  ■  -0.2  mm  the  variation  is  more  nonlinear,  especially 


the  crack  cli^M^g 
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near  y  -  0.  In  the  following,  (P^,  -  Pjj)  and  M  were  applied  at 
X  ■■  —0.76  mm. 

The  unit  load  solutions  F  ,  F  ,  d  ,  and  d  are 

X  y  X  y 

F^  »  9.36  X  10“^ 

X 

2  -1 

-  0.531  m  ^ 

X 

F^  -  0.0261 

y 

F^  -  -0.252  ra"^ 

-  1.40  X  10"^°  m-N~^ 

d^  -  1.17  X  10~®  N”^ 

y 

These  values  were  used  in  equation  5  for  any  cetobination  of  (P^  -  and  M 
to  obtain  and 

Differences  between  Gj.  and  Gj^  from  the  approximate  superposition 
analysis  and  the  geometrically  nonlinear  finite  element  analysis  can  be  traced 
mainly  to  two  sources:  (1)  nonlinear  variation  of  the  strains  through  the 
thickness  of  the  buckled  region  and  (2)  inaccuracy  in  determining  (P^  -  Pj^) 
and  M. 

By  using  (P^  -  Pjj)  and  M  from  the  geometrically  nonlinear  finite 
element  analysis,  the  effect  of  nonlinear  variation  of  the  strains  can  be 
exaitiined.  Figure  5  shows  that  this  effect  is  small. 

Figure  6  shows  chat  if  the  strength  of  materials  analysis  is  used  to  calcu¬ 
late  (P^  -  Pjj)  and  M,  the  difference  is  urach  larger.  Hence,  most  of  the 
difference  between  the  two  analyses  is  due  to  inaccuracy  in  determining 
(Pc  -  Pjj)  and  M.  But  the  general  trends  for  the  G^  variation  with 
delamination  length  and  load  are  predicted  very  well.  In  figure  6  the  curves 
for  the  two  analyses  seem  to  differ  (approximately)  by  a  constant  scale  factor- 
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A  direct  test  of  the  approxljnate  analysis  fcr  predicting  trends  is  to  use 
it  to  coalesce  the  curves  in  figure  6  into  a  single  curve.  Equations  (2)  and 
(3)  show  chat  (P^  -  and  M  can  be  expressed  as  functions  of  6.  Hence, 

from  equations  (5),  Gj  and  Gjj  are  functions  of  5.  Equations  (1),  (2), 
(3),  and  (S;!  show  that  for  constant  6,  P^  varies  as  a~^  and  Gj  and  G^j 
vary  as  a“^.  Hence,  plotting  a^Gj  vs  a^j  should  coalesce  the  curves  for 
various  delamination  lengths.  Figure  7  shows  that  the  data  for  five  delami- 
nat^on  lengths  (including  those  in  fig.  6)  do  coalesce  into  a  narrow  band 
around  a  single  curve.  Since  the  peak  values  of  Gj  for  various  lengths 
differ  by  more  than  two  orders  of  magnitude,  the  closeness  of  the  fit  suggests 
the  approximate  anaJ.ysis  is  accurate  for  predicting  trends.  Therefore,  all 
results  that  follow  are  obtained  with  the  approximate  superposition  analysis. 
Figure  7  also  shows  that  if  nonlinear  finite  element  results  are  available  for 
one  delamination  length,  the  values  for  other  lengths  can  be  estimated 
loonediately  • 

An  advantage  of  the  superposition  analysis  is  that  it  allows  a  problem  to 

be  dissected.  In  particular,  one  can  determine  the  relative  Importance  of  the 

loads  (P^  -  Pjj)  and  M  on  Gj  and  Gjj.  Figures  8  and  9  show  Gj  and  Gjj 

calculated  by  using  M  alone  and  by  using  (P^  -  and  M  in  combination. 

Although  intuition  might  suggest  that  only  the  peeling  action  caused  by  the 

moment  M  has  a  significant  effect  on  Gj,  figure  8  shows  that  (P_  -  P_) 

contributions  cannot  be  ignored.  Figure  9  shows  that  both  (P  -  P  )  and  M 

w  D 

are  also  important  when  calculating  Gjj, 

Parametric  Study 

The  effects  of  several  parameters  on  and  G^^  were  examined  using 
the  approximate  superposition  analysis.  The  parameters  were  initial  waviness, 
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delamination  length,  applied  load,  and  the  ratio  of  axial  to  bending  stiffness 
for  the  buckled  region. 

Initial  imperfections  in  the  form  of  simple  sinusoidal  wavinese  were 
assumed  (eq.  (6)) 


v(x) 

INITIAL 


-cos 


(6) 


where  v(x)  -  Distortion  in  the  y  direction.  When  a  column  is  initially 
wavy,  bifurcation  buckling  does  not  occur.  As  soon  as  load  is  applied,  the 
column  begins  to  deflect  laterally,  which  causes  interlaminar  stresses.  Hence, 
Gj  and  are  nonzero  as  soon  as  load  is  applied.  If  6^  "  0,  Gj  and  G^j 

are  zero  until  buckling  occurs.  However,  figure  10  shows  that  the  peak  value 
of  Gj  is  significantly  reduced,  even  for  very  small  imperfections.  In 
contrast,  figure  11  shows  that  G^^  is  hardly  affected  by  initial  waviness. 

Figures  6,  12,  and  13  show  the  effect  of  delamination  length  on  G^  and 
Gif.  The  shorter  delaminations  have  the  larger  values  of  peak  G^.  (i*e.,  Gj) 
(fig.  6).  However,  for  the  longer  delaminatlons  G^  becomes  nonzero  at  lower 
loads.  Figure  12  shows  that  after  only  a  little  delamination  growth,  G^ 
reaches  a  peak  and  decreases  rapidly  with  further  growth.  At  2a  ■  40-50  mm, 
the  crack  tip  closes  in  the  normal  direction  and  G^  is  identically  zero. 
Further  delamination  growth  causes  compressive  normal  stresses  to  develop  at 
the  crack  tip.  In  contrast,  Initially  increases  then  decreases  only 

slightly  to  a  constant  value  with  increased  delamination  length  (fig.  13). 

Note  that  G^^  is  typically  much  larger  than  Gj. 

Figures  6  and  14  illustrate  the  effects  of  applied  load  on  G^  and 
Gjj,  respectively.  The  mode  I  strain  energy  release  rate  G^  first 
increases  to  a  peak  value  (G^),  then  decreases  with  increasing  load 
(fig.  6).  In  contrast,  G^^  monotcr.lcally  increases  with  increasing 
load.  As  a  result,  and  G^^  do  not  usually  reach  peak  values 
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•t  the  ...e  derln*  .  fatigue  loed  cycle.  Futtheneote.  the  le«.  at  uhlch 
Cj  la  »e»l»u.  deeteaee.  with  lucreaalng  del«.ln.tlen  length. 

Since  Gj  Is  the  pe.eigia  value  of  0^  lor  a  given  delasduatlon 

length.  It  la  of  Intereat  ho.  .atles  with  delawlnatlnn  length.  The  first 
step  in  detenslnlng  thl.  variation  la  to  deteroUne  the  lateral  deflection  i 
corresponding  te  G^.  The  lateral  d.flectlon  J  la  oht.lned  hy  aolvlng 


3G, 


36 


(7) 


Eguatlon.  p,.  (5,,  .„g  are  eo.d.lned  to  obtain  the  governing 

equation  (eq.  (8)) 


<5(1  +  5^) 


D6 


1  -  S. 


Sa  + 


Sd 


(6  +  5.)  2 


+  -f- 


(8) 


Eguetlon  (8)  la  solved  Iteratively  for  S.  Once  J  ,a  detennlned.  can  be 
calculated  fro.  eduetlon  ,1).  (3,.  and  (5,.  gore  that  ?  la  mgependent 

of  delenlnatlon  length.  Earlier  It  was  .how.  that  for  eonatant  d.  G  and 


Gjj  vary  as  a"-^  and  P  varies  as  a'^.  Hence  G  v.rfd..  -4 

1  varies  as  a  and 


the  corresponding  applied  load  vanes  a.  a'^.  The  corresponding  value  of  0 

also  varies  a.  a^*.  These  ohservetlon,  will  he  of  .p«:lal  interest  later  when' 
examining  the  fatigue  data. 

The  last  para«ter  to  he  exanned  la  the  ratio  of  fle,urel-to-axUl 
stlffnee.  (l.c.,  D/S^)  of  the  delandnatcd  region.  The  buckling  load  for  the 
region  1.  linearly  related  to  B.  Prior  to  buckling,  the  load  In  the  dela.1- 
nated  region  1.  llne.rly  related  to  s„.  genre,  tor  thick  specl^n,  with  a 
thin  del«.lnated  region  the  applied  loed  which  cause,  buckllug  depends  on  the 
ratio  D/s„.  Delanlneted  regions  are  less  prone  to  buckle  It  they  possess  low 
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axial  stiffness  and  high  flexural  stiffness.  But  for  homogeneous  materials  or 
unidirectional  orthotropic  laminates,  this  ratio  is  simply 


(9) 


Hence,  the  applied  load  which  causes  buckling  is  independent  of  the  material 
properties  cf  the  buckled  region.  However,  for  multi-directional  laminates 
the  ratio  D/S^^  depends  on  both  the  lamina  properties  and  the  stacking 
sequence.  For  example,  the  value  of  for  a  [O^]  laminates  is  approximately 

1.4  times  that  for  a  [0/90/0]  laminate,  but  the  value  of  D  is  essentially 
the  same  for  both  laminates.  Consequently,  for  a  thick  laminate  with  a 
delamination,  buckling  occurs  at  a  lower  applied  load  if  the  delaminated  region 
is  rather  than  [O/90/O]. 


Comparison  of  Analysis  and  Fatigue  Data 

The  roles  of  and  in  delamination  growth  were  investigated  by 

comparing  calculated  values  of  G^  and  with  measured  growth  rates. 

Fatigue  data  for  six  specimens  from  reference  3  were  used  for  comparison  with 

analytical  results.  Three  of  the  specimens  were  tested  at  (P.p)  «  33  KN 

^  max 

and  another  three  at  (P«)  =  30  KN.  The  results  are  presented  in  figure  15. 

^  max 

^  decreased  rapidly  with  delamination  growth.  Both  curves  are  approximately 
dN 

linear  with  a  slope  of  -4,  hence 


/■ 


%  -  Ca"^  (10) 

oN 

da 

Figure  16  shows  —  vs.  the  maximum  values  of  Gj  and  G^^  during 
fatigue  cycling  at  a  maximum  compressive  load  of  33  KN.  Note  that  the  growth 
rate  is  largest  when  is  relatively  large.  But  slow  growth  is  observed 
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even  when  is  very  small.  The  mode  II  strain  energy  release  rate 

Gjj.  changes  little  after  Initial  delamlnatlon  growth.  Since  G 
remains  large,  delamlnatlon  growth  likely  is  driven  by  G^^. 

Two  delamlnatlon  growth  criteria  were  examined  to  determine  whether  they 

could  predict  the  observed  growth  rates.  The  first  growth  criterion  examined  is 
given  by 


where  Z  and  n  are  constants.  Eqn.  (11)  m  evaluated  at  the  point  In  the 

load  cycle  when  G^  is  maximum.  The  load  range  and  delamlnatlon  lengths  are 
such^that  for  almost  the  entire  test,  the  maximum  G^  during  each  load  cycle 
is  Gj,  which  was  obtained  by  solving  equation  (7).  Recognizing  that  G 
decreases  as  a'^  as  shown  earlier,  equations  (10)  and  (11)  can  be  solved  for 
n;  the  result  is  n  .  1.  This  is  in  strong  contrast  to  published  values  of 
n  S  15-20  for  double  cantilever  beam  fatigue  tests  (ref.  3).  Apparently, 
equation  (11)  is  not  a  valid  growth  criterion  for  the  specimens  considered. 

Next  a  growth  criterion  was  considered  which  includes  both  G^  and  G 
If  we  assume  there  is  no  synergistic  interaction  of  G^  and  G^^  (i.e.,  the 

effects  are  separable) ,  then 

(12) 

where  f^  and  f^  are  functions  of  G^  and  G^^  respectively. 

From  the  double  cantilever  beam  data  just  discussed,  we  know  that  f, 

1.  .xtre..!,  to  (i...,  ^  -  cj’).  since  0,  Oecreese.  ' 

rapidly  with  increasing  "a”,  f^  must  also  decrease  extremely  fast  as  "a" 
Increases.  In  fact,  f^  would  not  contribute  noticeably  to  da/dN  after 
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the  initial  groirth.  Bence  delaainatlon  growth  appears  to  be  driven  by 
alone.  Accordingly,  it  was  assuaed  that  the  growth  criterion  should  be 
evaluated  when  la  maximum,  l.e.,  at  peak  load.  However,  earlier  It 

was  shown  that  for  long  delamlnatlona  the  crack  tip  closes  and  produces  com¬ 
pressive  Oy  stresses  when  the  cyclic  load  Is  maximum.  The  compressive 
stress  probably  reduces  the  effect  of  G^j  on  delamlnatlon  growth,  but  It 
was  not  clear  how  to  account  for  this  stress.  Two  approaches  were  tried: 

(1)  Ignore  the  compressive  uormal  stress  and  set  f ••  0  when  the  tip 
closes  or  (2)  let  f^  take  on  negative  values  after  the  crack  tip  closes. 

If  we  choose  to  set  “  0  when  the  crack  tip  closes,  then 

^  “  ^2^‘^II>  (13) 

for  virtually  the  entire  test.  Figure  14  showed  that  0^.^.  first  Increased 
then  decreased  slightly  as  the  delamlnatlon  extends.  In  the  experiments  two 
load  levels  were  used:  -  33  KN  and  30  KN.  Figure  14  shows  that  for 

2a  >  25  mm,  the  minimum  value  of  G  for  (p^)  -  33  KN  Is  greater  than 

the  maximum  value  of  G^^  for  “  30  KN.  Hence,  equation  (13)  would 

predict  that  for  2a  >  25  m,  the  minimum  da/dN  for  the  higher  load  should 
exceed  the  maximum  da/dN  at  the  lower  load.  Figure  16  shows  this  Is  not  the 
case.  Hence,  equation  (13)  Is  not  valid. 

If  we  select  a  function  fj^  that  becomes  negative  when  the  crack  tip 
closes,  then  we  (analytically)  allow  compressive  normal  stresses  at  the  crack 
tip  to  retard  delamlnatlon  growth  due  to  G^j.  Since  the  compressive  crack  tip 
stresses  Increase  as  the  delamlnatlon  grows,  such  a  function  f^^  would  predict 
a  decrease  In  growth  rate  with  Increased  delamlnatlon  length.  Although  this 
prediction  agrees  with  the  data  trend  In  figure  15,  more  tests  are  needed  to 
verify  or  disprove  this  Interpretation. 
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Despite  the  coiaplexlty  of  the  grovth  behavior,  two  trends  were  clear: 

(1)  high  growth  rates  were  observed  only  when  ras  large,  and  (2)  slow 
growth  was  observed  even  when  Gj  was  negligibly  small;  apparently,  Gjj 
alone  can  drive  delanlnation  growth. 

CONCLUSIONS 

Analysis  and  experiments  were  used  to  study  instability-related  delamina- 
tion  growth  in  a  fatigue  specimen  with  a  through-width  delamination.  To  per¬ 
form  the  analysis  an  approximate  superposition  analysis  was  developed.  The 
analysis  expresses  Gj  and  ®II  in  closed  formg  which  can  be  used  easily  to 
determine  the  effects  of  various  parameters*  The  analysis  agreed  very  well 
with  more  rigorous  solutions. 

The  response  of  the  delaminated  laminate  to  applied  loads  was  found  to  be 
very  complex*  Key  observations  are  listed  below* 

(1)  Gjj  is  generally  much  larger  than  Gj* 

(2)  Gj  and  Gjj-  usually  reach  their  peak  magnitudes  at  different  points 

in  a  fatigue  cycle*  Gjj  always  reaches  its  peak  value  at  maximum 
load* 

(3)  High  delamination  growth  rates  were  accompanied  by  large  values  of  Gj 

(4)  Slow  growth  rates  were  observed  even  when  Gj  was  negligibly  small. 

This  growth  apparently  was  due  to  a  large  value  of 
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APPENDIX 


A  strength  of  materials  analysis  is  described  herein  for  the  configuration 
in  figure  1. 

The  configuration  is  subdivided  into  four  regions »  as  shown  in  figure  2. 
Because  of  symzaetry,  only  half  of  the  laminate  is  modeled*  The  laminate  is  of 
width  b*  The  following  assumptions  are  made: 

1.  Regions  B  and  C  are  perfectly  bonded*  Regions  A  and  D  are  unbonded* 

2*  Regions  A»  B»  and  C  have  constant  axial  strain.  Hence»  the  force* 
displacement  relations  are  those  for  a  simple  rod  subjected  to 
axial  load* 

3*  Region  D  has  zero  slope  at  both  ends* 

A.  Region  D  has  an  initial  sinusoidal  imperfection  of  peak  magnitude 
5^*  The  initial  shape  is  given  by 

where  v(x)  •  the  distortion  in  the  y  direction. 

To  describe  the  nonlinear  behavior  of  region  D,  equations  (A2)  and  (A3)  for 
post-buckling  of  a  column  were  used* 


P  6 

D  2  6  +  (S 
a  o 


(ref.  4)  (A2) 


(ref.  5) 


(A3) 


where  5,  a,  a,  and  are  peak  lateral  deflection,  axial  length  before  and 
after  deformation,  and  load,  respectively.  Equations  (A2)  and  (A3)  were 
derived  using  strength  of  materials  analysis  of  a  column. 
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To  combine  regions  A,  B,  p,  and  D,  equilibrium  and  compatibility 

conditions  must  be  considered.  The  equUibrium  condition  for  the  axial  force 
is 


■•a  +  ‘■o  ■  +  fc  ■  '■t 


(A4) 


Compatibility  requires  the  shortening  of  regions  A  and  D  to  be  identical. 
Hence, 


a  -  a 


(A5) 


Equations  (A2)  to  (A5)  can  be  combined  to  obtain  the  governing  equation  for 
the  laminate • 


TT  S. 


r  1  2  X  1 

(5^  +  266  )  +  - P  . 
[16  o^  6  +  6^J 


(A6) 


For  a  specified  load  P^,  equation  (A6)  is  solved  using  a  Newton-Raphson 
technique  to  obtain  6,  can  then  be  calculated  using  equation  (A2).  From 

static  equilibrium,  the  moment  acting  on  the  delaminated  region  at  the  crack 
tip  is 

The  force  Is  found  from  rule  of  mixtures  as 

p  .-i£_  p 

c  (A8) 
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Fig.  1—  Local  buckling  of  laminate  with  through-^dth  delamination. 
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OEUMINATION 


Compaction  of  results  from  geometrically  nonlinear  analysis.  (Curve 
Is  visual  fit  of  data.) 


0  mm 


0  mm 


ORIGINAL  PAGE  13 
OF  POOR  QUALITY 


J _ I _ L_ _ J _ I _ I 


APPLIED  COMPRESSIVE  LOAD. 


36  kN 


ORIGINAL  PAGE  IJ 
OF  POOR  QUALITY 


m  % 


•  *  ^ 


I  cr\l  O 

I  cr\i  cr\ 


f 


CM 


-.8 


o 

M 

M 

oo 

E 

E 

«h 

O 

g 

<K3 

42 

CM 

4J 

00 

p 

— 

rr 

Q> 

25  mm 


35 


Fig.  14-  Effect  of  applied  compressive  load  on  Gjj  (2a  -  25  mm). 


